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ABSTRACT 

We present an online catalog of distance determinations for 4781 K giants, 
most of which are members of the Milky Way's stellar halo. Their spectra from 
SDSS/SEGUE provide metallicities with accuracies A[Fe/H] ^ ±0.2 dex and 
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giant-dwarf distinction. The distance moduli are derived from a comparison of 
each star's apparent magnitude with the absolute magnitude of empirically cali- 
brated color-luminosity fiducials, at the observed {g — r)^ color and spectroscopic 
[Fe/H]. We employ a probabilistic approach that makes it straightforward to 
propagate the errors in metallicities, magnitudes, and colors properly into dis- 
tance uncertainties. We also fold in prior information about the giant-branch 
luminosity function and different metallicity distributions of the SEGUE K-giant 
targeting sub-categories. We show that the metallicity prior plays little role in 
the distance estimates, but that neglecting the luminosity prior would lead to a 
systematic distance modulus bias of up to 0.2 mag. We find a median distance 
precision of 12%, with distance estimates most precise for the least metal-poor 
stars near the tip of the red-giant branch. We use globular and open clusters 
to verify the precision and accuracy of our distance estimates. The stars in our 
pubhcly available catalog are up to 110 kpc distant from the Galactic center, 
with 270 stars beyond 50 kpc, forming the largest sample of distant tracers in 
the Galactic halo. 

Subject headings: galaxies: individual(Milky Way) - Galaxy: halo - stars: K 
giants - stars: distance 



1. Introduction 



K giants have lorig been used t o map the Milky Way's s tellar halo ( lBondlll980l : lRatnatunga fc Bahcall 



19851 : iMorrison et al.lll990l l2000l : IStarkenburg et al.ll2009l ). In contrast to blue horizontal- 
branch (BHB) and RR Lyrae stars, giant stars are found in predictable numbers in old 
populations of all metallicities, and at the low metallicities expected for the MW halo they 
are predominantly K giants. At the same time, their high luminosities (M^ ~ 1 to —3) make 
it fea sible to study th em with current wide-field spectroscopic surveys to distances of > 100 
kpc flBattaglial l2007l) . The Sloan Extension for Galactic Understanding and Exploration 
(lYanny et al.ll2009l . SEGUE), which has entered the second stage (SEGUE-2, Rockosi et al. 
in preparation) , specifically targeted K giants for spectroscopy. For simplicity, henceforth we 
refer to SEGUE and SEGUE-2 collectively as simply SEGUE. The SEGUE data products 
include sky positions, radial velocities, apparent magnitudes, and atmospheric parameters 
(metallicities, effective temperatures, and surface gravities), but no distances. 

Distance estimates to kinematic tracers, such as the K giants, are i ndispensable for stud- 



ies of Milky Wa y halo dynamics - such as estimates of the halo mass (iBattaglia et al.l 12005 



Xue et al.ll2008l ). for exploring the formation of our Milky Way (e.g., probing velocity-position 
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correlations, IStarkenburg et al.ll2009l : IXue et al.ll201l[ ). and for deriving the metallicity pro- 
file of the Milky Way stellar halo. All of these studies require not only unbiased distance 
estimates, but also a good understanding of the distance errors. However, unlike 'standard 
candles' (i.e., BHB and RR Lyrae stars), the (intrinsic) luminosities of K giants vary by two 
orders of magnitude, with color and luminosity depending on stellar age and metallicity. Fur- 
thermore, the lower (less luminous) part of the giant branch is 'steep' in the color-magnitude 
diagram (CMD), particularly for stars of low metallicity. For this reason, it turns out to be 
more difficult to obtain precise distances for K giants than for, say, main-sequence stars. 

The most immediate approach for estimating a distance to a K giant with color c and 
metallicity [Fe/H] ( e.g., from SDSS/SEGUE) is to simply look up its expected absolute 
magnitude M in a set of observed cluster fiducials, M(c, [Fe/H]). Comparison with the 
apparent magnitude then yields the distance modulus (denoted by VAi) and distance. In 
practice, this simple approach has two problems. First, care is required to propagate the 
errors in metallicities, magnitudes, and colors properly into distance uncertainties. Secondly, 
such an approach does not immediately incorporate external prior information, such as the 
luminosity function along the red giant-branch (RGB) and the overall metallicity distribution 
of the stellar giant population under consideration. Because the luminosity function along 
the RGB is steep, a nd there are many more low er luminosity stars than high luminosity 



stars (n(L) ~ L ISandquist et al.lll996l Il999[ ). an estimate of the absolute magnitude 



M(c, [Fe/H]), is more likely to produce an overestimate of the luminosity, and therefore an 
overestimate of the distance. Analogously, there are few extremely metal-poor or metal- 
rich stars observed in the halo; hence, a very low estimate of [Fe/H] is more likely an 
underestimate of the metallicity of an (intrinsically) less metal-poor star, and so the estimate 
of the absolute magnitude will lead to an overestimate the luminosity, and therefore an 
overestimate of the distance. Thus, to exploit K giants such as those from SDSS/SEGUE for 
various dynamical analyses, an optimal way to determine an unbiased distance probability 
distribution for each sample star is crucial. 

A general probabilistic framework to make inferences about parameters of interest (e.g., 
distance moduli) in light of direct observational data and broader prior information is well- 
established. It has been applied in a wide variety of circumstances and recently also ap- 
plied to the distance det e rminations for stars , i ncluding giant star s in th e RAVE survey 
teurnett &: BinnevI lioiol : iBurnett et all boilj ). iBurnett fc BinnevI J2OI0I ) described how 
probability distributions for all the 'intrinsic' parameters (e.g., true initial metallicity, age, 
initial mass, distance, and sky position) can be inferred according to Bayes' theorem, draw- 
ing on the star's observables and the errors thereon. Here we focus on a somewhat more 
restricted problem: the distances to stars on the RGB, which we can presume to be 'old' 
(> 5 Gyr). Like any Bayesian approach, our approach is optimal in the sense that it aims 
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to account for all pertinent information, can straightforwardly propagate the errors of the 
observables to distance uncertainties, and should avoid systematic biases in distance esti- 
mates. This approach also provides a natural framework to account for the fact that distance 
estimates will be less precise for stars that fall onto a 'steep' part of the color-magnitude 
fiducial, presumably metal-poor stars on the lower portion of the RGB. 

The goal of this paper is to outline and implement such a Bayesian approach for esti- 
mating the best and unbiased probability distribution of the distance moduli VAi, for each 
star in a sample of 4781 K giants from SDSS/SEGUE. This distribution can then be char- 
acterized by the most probable distance modulus, PA^peak, and the central 68% interval, 
AVAi. At the same time, this approach of course also yields estimates for the absolute 
magnitude M, heliocentric distance d, Galactocentric distance rcc, and their corresponding 
errors. 

In the next section, we introduce the selection of the SEGUE K giants and their ob- 
servable properties. In §3, we describe a straightforward (Bayesian) method to determine 
the distances. The results and tests are presented in §4. Finally, §5 presents our conclusions 
and a summary of the results. 



2. Data 



SDSS and its extensions use a dedicated 2.5m tele s cope (IGunn et aLll2006h to obtain 



uqriz 



2002: Pier et al. 



imaging (iFukugita et al. 



2003 



1996 



Eisenstein et al. 



2011 



Gunn et all Il998l : I York et all I2OOOI : IStoughton et al. 



Smee, S. A., et al. 2012, AJ, submitted) and 
resolution (defined as R = A/AA) ~1800 spectra for 640 objects in a 7 sq deg field simulta- 
neously. SEGUE (Sloan Extension for Galactic Understanding and Exploration), one of the 
key projects of SDSS-II and SDSS-III, obtained some 360,000 spectra of stars i n the Galaxy, 



selected to explore the nature of stellar populations from 0.5 kpc to 100 kpc (lYanny et al. 



2OO9I . and Rockosi et al in preparation). Data from SEGUE is a significant part of the ninth 
public data release, DR9 (Ahn, C., et al. 2012, ApJS, submitted). 

SDSS DR9 delivers estimates of Teg, log g, and [F e/H] from an upd a ted and improved 
version of the SEGUE Stella r Parameter Pipeline fSSPP- lhee et al.ll2008al lbl: I AUende Prieto et al, 



20081 : ISmolinski et al.ll201ll ). The SSPP processes the wavelength- a nd flux-cahbrated spec 



tra generated by the standard SDSS spectroscopic reduction pipeline (jstoughton et al. 2002), 
obtains equivalent widths and/or line indices for more than 80 atomic or molecular absorp- 
tion lines, and estimates Teg, log g, and [Fe/H ] through the application of a number of 
complementary approaches (see lLee et al.ll2008al . for detailed description of these techniques 
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and Rockosi et al. in preparation for the changes and improvements of the SSPP). 

Targeting K giants has been an important part of SEGUE's scientific aims, particularly 
in SEGUE-2. Because of the intrinsically low fraction of K giants in the Galaxy, SEGUE 
adopted sophisticated criteria mainly based on ugriz photometry and proper motions to 
preselect likely K-giant candidates for fiber-fed spectroscopic observation. 

SEGUE targets K-giant candidates with [g — r)o between 0.5 and 1.3, {u — g)^ between 
0.5 and 3.5 (see Figured]), and proper motions smaller than 11 mas/year. The K giants are 
designed to be selected in three sub-categories, named l-coloi0 K giants, red K giants, and 
proper motion K giants. However, besides the three sub-categories, about 20% of our giant 
sample comes from stars which do not pass the target criteria described above (hereafter we 
referred to as serendipitous K giants.). The majority of serendipitous K giants were originally 
targeted as K giants, but as photometric reductions have improved, their ugr colors changed 
so that they are no longer in the target region. There are also a number of serendipitous 
K giants from different target-selection classes such as low-metallicity stars (which use a 
stronger cut on 1 color and no proper motion criteria) and G and K dwarfs (d esigned to 



study the disk, with only a simple cut in g — r). Figure 10 of lYanny et al.l (120091 ) shows the 
regions of the u — g/g — r plane occupied by K giants: each target category focuses on a 
particular region. The 1-color K-giant category uses the metallicity sensitivity of the u — g 
color for the bluer part of our color range to select metal-poor K giants. The other two 
categories focus on the redder stars with [g — r)o > 0.8. The red K-giant category selects 
those stars whose luminosity leads them to fall above the locus of foreground stars, while 
the proper motion K-giant region is where the K giants are found in the locus of foreground 
stars. Here only a proper motion selection can be used to cull out the nearby dwarf stars 
with apprecia ble proper moti o n. W e further limit the sample by requiring that the reddening 



estimate from ISchlegel et al.l ( 1l998l ) for each star E(B-V) is less than 0.25 mag. 



In addition to the log g estimate provided by the pipeline, we employ a spectroscopi c 



index in the region of the Mgb/MgH features, similar to that of [Morrison et al.l ( 120031 ). 
which is a strong indicator of luminosity. The SSPP continuum corrects the spectra using 
a polynomial, and this process can remove much of the strong MgH feature seen in dwarfs; 
thus we work with flux-corrected but not continuum-corrected data. We calibrate this index 
using spectra taken of globular and open cluster giants. This is described in full in Ma et al. 
(in preparation). 

Using colors, reddening, [Fe/H], and log g values from DR9, we find 14,025 field candi- 



^1 — color = — 0.436u + 1.129g — 0.119r — 0.574i + 0.1984 and i s a ph otometry metallicity indicator for 



stars in the color range 0.5 < [g — r) < 0.8. Reference: iLenz et al.l (|199 
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dates which satisfy our K-giant criteria and have good photometry (i.e., color errors are less 
than 0.04 mag). We describe a further culling of the sample in §3.4, where we also describe 
the fiducials used to obtain distances to eliminate stars that could either be on the RGB or 
in the red clump. The sample has a typical color error of ±0.02 mag, a typical [Fe/H] error 
of ±0.2 dex, and a typical log g error of ±0.5 dex. The error of [Fe/H] for each K giant used 
in this paper is calibrated using cluster data plus repeat observations, which depends on the 
signal-to- noise ratio of the spectrum, as described in detail in Ma et al. (m preparation). 



3. Probabilistic Framework for Distance Estimates 



Our goal is to obtain the posterior probability distribution (pdf) for the distance modulus 
of any particular K-giant star, after accounting for (i.e., marginalizing over) the observational 
uncertainties in apparent magnitudes, colors, and metallicities (m, c, [Fe/H], Am, Ac, A[Fe/H]), 
and after including available prior information about the K-giant luminosity function, metal- 
licity distribution, and, possibly, the radial halo-density profile. 



3.1. Distance Modulus Likelihoods 

We start by recalling Bayes theorem, cast in terms of the situation at hand: 

P({m, c, [Fe/H]} | VM) Pprior {VM) 



F{VM I {m, c, [Fe/H]}) 



P({m, c, [Fe/H]}) 



Here, P(T>Ai \ {m, c, [Fe/H]}) is the pdf of the distance moduli {VAi = m — M), and de- 
scribes the relative probability of different VJH, in light of the data, {m, c, [Fe/H]} (we use { 
} to denote the observational constraints, i.e., the estimates and the uncertainties for the ob- 
servable quantities). P({m, c, [Fe/H]} | DM) is the likelihood oi VM (e.g., ^{VM)), and 
tells us how probable the data {m, c, [Fe/H]} are if VM were true. The term Ppnor{VM) is 
the prior probability for the VM, which reflects independent information about this quantity, 
e.g., that the stellar number density in the Galactic halo follows a power law of r~^. The 
term P({m, c, [Fe/H]}) is a non-zero constant. 

So, the probability of the VM for given star is proportional to the product of the 
likelihood of VM and the prior probability of VM (e.g., Eq 2). 



P{VM I {m, c, [Fe/H]}) oc ^{VM) Pprior [VM] 



(2) 
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The prior probability for the DM in Eq. 2 can be incorporated independently, and the 
main task is to derive ^(T>A4) . In deriving ^(T>Ai) , we must in turn incorporate any prior 
information about other parameters that play a role, such as the giant-branch luminosity 
function, Pprior(M), or the metallicity distribution of halo giants, Ppi.ior([Fe/H]). This is done 
via 



^{VM) = J J p({m, c, [Fe/H]} | DTU, M, FeH)ppHor(M)ppHor(FeH)dM dFeH. (3) 

Here we use FeH to denote the metallicity of the model, while we use [Fe/H] for the 
observed metallicity of the star. 



3.2. Observables and Priors 

The direct observables we obtain from SEGUE and SSPP are the extinction corrected 
apparent magnitudes, colors, metallicities, and their corresponding errors (r^, {g — r)o, [Fe/H], Arg, 
A{g — r)o, A[Fe/H]). Figure [2] shows that the most common stars are the intrinsically fainter 
blue giants, as we would expect from the giant-branch luminosity function. Hereafter, we 
keep using c and m instead of {g — r)^ and tq for convenience and generality in the expression 
of the formulas. 

In our analysis, we can and should account for three pieces of prior (external) information 
or knowledge about the RGB population that go beyond the immediate measurement of the 
one object at hand: Pprior(^^A^), Pprior(M) and Pprior([Fe/H]) (shown as Eq. 1 and Eq. 3). 

The prior probability of V Ai reflects any infor matio n on the r a dial d ensity profile of 



the Milky Way's stellar halo. IVivas &: ZinnI (|2006|) and ISell et al.l (|2008|) indicated that 



the radial halo stellar density follows a power law p oc r", with the best value of a = 
—3, and reasonable values of —2 > a > —4; this implies Pprior(^^A^)d'DA^ = p47rr^dr, 
Pprior(^^A^) + ^ Quite fortuitously, the prior probability for VAi turns out to 

be fiat for the radial stellar density of a power law of p oc r~^. Given that the =Sf (PA^) 
approximatively follows a Gaussian distribution with the mean of PA^o ^'^^ the standard 
deviation of AVAi (here AVAi is the error of VAi), the Pprior(^Al) will shift the estimate 
ofVMo by i^±^^2i£i°(Ar'Al)2, but with basically no change in AVM. Therefore, the shifts 
in the mean VAi caused by Pprior('^'A^) can be neglected for values of — 2 > a > — 4 (i.e. 
^2+^^f^{AVMf < AVM). 

The prior probability of the absolute magnitude M can be inferred from the nearly 
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universal luminosity function of the giant branch of old stellar population s . Spe cifically, we 



derive it from the globular clusters M5 an d M30 (ISandauist et a. 



the Basti theoretical luminosity function (jPietrinferni et al.l |2004| ). Figure [3] (top panel) 



19961 . 119991 ). and from 



shows that the luminosity functions for the RGBs derived from the globular clusters in the 
different bands are consistent with each other, and also consistent with the Basti theoretic 
luminosity functions for the metal-rich and metal-poor cases. All the luminosity functions 
follow linear functions with similar slope, k = 0.32, as a result that the luminosity functions 
are insensitive to changes in the metallicity and wave bands. According to p(M)dM = p(L)dL 
and M ~ — 2.51ogL, the luminosity function p(M) oc lO'^'^ means p(L) oc L"^'^'^"^. We 
conclude that the luminosity function for the giant branch follows p(L) oc L~^'^, shown as 
the dashed line in Figure [31 

Our prior probability for [Fe/H] results from an empirical approach. In the SEGUE 
target selection, the K giants were split into four sub-categories: red K giants, 1-color K 
giants, proper motion K giants, and serendipitous K giants. This suggests that one should 
adopt the overall metallicity distribution of each sub-category as the [Fe/H] prior for any one 
star in this sub-category (Figure Hj). Figure H] shows the metallicity distribution along with 
apparent magnitude in the upper panel, and the four [Fe/H] priors in the lower panel. It can 
be seen that the four sub-categories have different metallicity distributions. For a star that 
has approximately the mean metallicity, this prior should leave ^(VAi) unchanged, because 
the individual metallicity error is smaller than the spread of the pprior([Fe/H]). However, 
for a star of seemingly very low metallicity, the prior implies that this has been more likely 
a metallicity-underestimate of a (intrinsically) less metal-poor star, which would mean an 
overestimated distance modulus. 



3.3. Color-Magnitude Fiducials 



To obtain distance estimates, we determine an estimate of the absolute magnitude 
of each star using its {g — r)^ color and a set of giant-branch fiducials for clusters with 
metallicities ranging from —2.38 to -1-0.39, an d then use the star 's apparent magnitude 
(corrected for reddening using the estimates of ISchlegel et al.lll998l ) to obtain its distance. 
We prefer to use Gducials, rather than model-isochrone giant branches wherever possible, 
because isochrone giant branches cannot reproduce cluster fiducials with sufficient accuracy. 

We derived such fiducials, using the globular clusters M92, M13, and M71, and the open 
cl uster NGC 6791 , which have accurate ugriz photometry from the DAOphot reductions 
of lAn et al.l (j2008[ ). For stars brighter than g ~ 14.5, the SDSS imager saturates, so we 
supplemented the giant-branch fiducials for the three globular clusters using the u'g'r'i'z' 
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m 



Tucker et al. 



ph otometry of IClem et al.l ( 120081 ) . We transformed to ugriz using the transformations given 



( 120061 ) ■ We give our adopted values of [Fe/H], reddening, and distance 
modulus for each cluster in Table 1. In addition, we changed the locus of the M71 fiducial 
slightly by only using stars which were known members. We supplemented the fiducials with 
a sol ar metallicity giant branch from the Basti a-enhanced isochrones ( iPietrinferni et al. 
20041 ). Figure [5] shows the four fiducials and one theoretical isochrone. The color at a given 



M and [Fe/H], c(M, [Fe/H]) can then be interpolated from these color-magnitude fiducials. 

Given the sparse sampling oi M - [g — r)o space by the four isochrones, we need to 
construct interpolated fiducials. We do this by quadratic interpolation of c(M, [Fe/H]), 
based on the three nearest fiducial points in color, and construct a dense color table for 
given M and [Fe/H], which will be used for Eq. 4. Extrapolation beyond the metal-poor 
and metal-rich boundaries and the tip of the RGB would be poorly constrained. Therefore, 
we use these limiting fiducials instead for the rare cases of stars with [Fe/H] < —2.38 or 
[Fe/H] > +0.39. Table 3 shows an interpolated fiducial with [Fe/H] = —1.18; the entire 
catalog of 20 interpolated fiducials with [Fe/H] ranging from —2.38 to +0.39 is provided 
in the electronic edition of the Journal. While there is more than one way to interpolate 
the colors, such as quadratic or piecewise linear, we have checked and found that different 
interpolation schemes lead to an uncertainty less than ~ ±0.02 mag in color, due to the 
sparsity of the fiducials. 



3.4. Red Giant-Branch Stars vs. Red-Clump Giants 

In addition, we have chosen not to assign distances to stars which lie on the giant 
branch below the level of the horizontal branch (HB). This is because the red horizontal- 
branch or red-clump (RC) giants in a cluster have the same color as these stars, but quite 
different absolute magnitude, and the SSPP log g estimate is not sufficiently accurate to 
discriminate between the two options. We derive a relation between [Fe/H] and the {g — r)o 
color of the giant bran ch at the level of the horizontal branch using eight clusters with 



ugriz photometry from I An et al.l (120081 ). with cluster data given in Table 2. The [Fe/H] 
and {g — r)^^ for the HB/RC of the clusters follow a quadratic polynomial, {g — r)^^ = 
0.086[Fe/H]^ + 0.38[Fe/H] + 0.96, as shown in Figure [2j We then use this polynomial and its 
[Fe/H] estimate to work out, for each star, whether it is on the giant branch above the level 
of the HB; it turns out that more than half of the candidate K giants fall into the region of 
RGB - HB ambiguity. 



Ultimately, we are left with 4781 stars that appear to be clearly on the RGB, above the 
level of the HB, according to the SSPP log g and the Mg index check. Of these, 195 satisfy 
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the target criteria for red K giants, 463 the criteria for proper motion K giants, and 3231 
the 1-color K-giant criteria. Another 892 were serendipitous identifications: stars targeted 
in other categories which nevertheless were giants. Figure |2] shows the distribution of the 
apparent magnitudes, ro, and metalhcities, [Fe/H], along with the color, [g — r)o. 

Besides the contamination from HB/RC stars, we need to consider possible contamina- 
tion of our sample by asymptotic giant-branch (AGB) stars, because it is not possible for 
us to distinguish between RGB and AGB stars with our spectra. While the difference in 
absolute magnitude can be large (varying from ~ 0.8 mag, implying a 40% distance underes- 
timate, at the blue end of our giant color range to a negligible difference at the giant-branch 
tip), the proportion of our giants which are on th e AGB is relatively small. We used both the 
luminosity function of ISandquist fc Boltd (j2004( ) for the globular cluster M5 and evolution- 
ary tracks from Basti isochrones for old populations of metallicity -2.6 and -1.0 to estimate 
the percentage of stars which are on the AGB. We find that for the most metal-poor stars, 
around 10% will be AGB stars, while for stars with [Fe/H] close to -1.0 the fraction is near 
20%. 



3.5. Implementation 

For any given star, the observables are its apparent magnitude and associated error, 
(mi,Ami), its color and error (q, Aq), and its metallicity and error ([Fe/H]j, A[Fe/H]j) . 
The VAi and the data are linked through the absolute magnitude M via: mi = M -|- VA4i 
and the fiducial c(M, FeH), which we presume to be a relation of negligible scatter. Now we 
can incorporate the errors of the data and the specific priors on the stellar luminosity and 
metallicity distribution when calculating J^'(T>Ai^ (see Eq. 3). 

In practice, the errors on color, apparent magnitude, and metallicity can be approxi- 
mated as Gaussian functions, in which case p({m, c, [Fe/H]} | DA^,M, FeH) (see Eq. 3) is 
modeled as a product of Gaussian distributions with mean and Delta (q, Aci), (rrij. Ami), 
and ([Fe/H].,A[Fe/H].): 



p({m, c, [Fe/H]}i | VM, M, FeH) 



1 



'271 A rrii 



exp 



{VM + M 



2(Am,) 



27rAci 
1 



exp 



X 



27rA[Fe/H] 



exp 



:c(M,FeH)-Ci)n 
2(Aci)2 / 
{FeH-[Fe/R]f - 
2(A[Fe/H]J2 ^ 



(4) 



For Equation 3, we use the priors Ppnor (M) , based on the luminosity function of the giant 
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branch, p(L) oc L (Figure [3]), and Pprior([Fe/H]) , based on the metalhcity distributions 
of the K-giant sub-categories (Figure Hj). 

For any K giant with {n\, [Fe/H].}, we can then calculate ^(T>A4^ by computing 
the integral of a bivariate function (Eq. 3) over dM and dFeH, using iterated Gaussian 
quadrature. The best estimate and the error in T>A4 are then estimated using the peak of 
the likelihood of VM (PA^peak) and its central 68% interval ( P^84%-pa^i6% 

To speed up the determination of the integral in Eq. 3, we look up c(M, [Fe/H]) in a 
pre-calculated, finely sampled color table instead of an actual interpolation. This approach 
can provide a consistent c for given M and [Fe/H], if the pre-prepared color table is suitable. 
We use a color table, c(M, [Fe/H]), of size 7000 x 8000, with -3.5 < M < 3.5 and -3.58 < 
[Fe/H] < +0.56. 



4. Results 

4.1. Distances for the SDSS/SEGUE K giants 

From §3.4, we obtain the Bayesian estimates of the distance moduli and their uncer- 
tainties from the peak of ^(VAi) and its central 68% interval. At the same time, we 
obtain estimates for the intrinsic luminosities by = rg — PA^peak, distances to the Sun by 
d = 10 5 , and Galactocentric distances by taking the Sun at = 8.0 kpc. This results 
in the main entries in our public catalog: the best estimates of the distance moduli and their 
uncertainties {VMpeak, ADM), distance moduh at (5%, 16%, 50%, 84%, 95%) confidence of 
J!f{T>M) {'DM5%,'DMi6%,'DM5o%,T^Msi%, T>Mq5%), distances to the Sun, and their errors 
(d. Ad), Galactocentric distances and corresponding errors {tqc, Arcc), and the absolute 
magnitudes along with the errors (M^, AM^) for 4781 K giants. In addition, we also include 
the sky position (RA, Dec), extinction-corrected apparent magnitudes, tq, colors, {g — r)Q, 
heliocentric radial velocities, RV, and corresponding errors, as well as the SSPP atmospheric 
parameters, Tefr, [Fe/H], log g, in the catalog for convenient use. Table 4 shows an example 
of the on-line table of the K giants. The complete version of this table is provided in the 
electronic edition of the Journal. 

Figure [7] illustrates the overall properties of the ensemble of distance estimates. The top 
two panels show the mean precisions of 12% in Ad/d and ±0.25 mag in these panels 

also show that the fractional distances are less precise for more nearby stars, because they 
tend to be stars on the lower part of the giant branch, which is steep in the color-magnitude 
diagram, particularly at low metallicities. 
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The bottom panel of Figure [7] shows that the mean error in is ±0.25 mag, and that 
faint giants have less precise intrinsic luminosity estimates. Figure [8] (upper panel) shows 
the distribution of K giants on the CMD. There are more stars in the lower part of the giant 
branch, which is consistent with the prediction of the luminosity function. The lower panel 
of Figure [8] shows that stars in the upper part of the RGB have more precise distances than 
those in the lower part of the CMD, because the fiducials are much steeper near the sub-giant 
branch. This is equivalent to the fact that the fractional distance precision is higher for the 
largest distances (see Figure [7]). 

The giants in our sample stars lie in the region of 5 — 110 kpc from the Galactic center 
and 5 kpc from the disk. There are 1478 stars beyond 30 kpc, 27 stars beyond 5 kpc , 
and 43 stars beyond 80 kpc (cf. The 5 red giants be yond 50 kpc in iBattaglia et al.l ( 120051 ) . 
16 halo st a rs beyond 8 kpc in iDeason et al.l (120121 ) and no BHB stars beyond 80 kpc in 
Xue et al.l (120081 . 1201 ll )). Our sample comprises the largest sample of distant stellar halo 
stars with measured radial velocities and distances to date. 



4.2. The Impact of Priors 

In this section we briefly analyze how important the priors actually were in deriving the 
distance estimates. For each star, the evaluation of Eq. 3 using Eq. 4 and the interpolated 
fiducials results in ^i(DM^ (Figure [6]), i.e., the likelihood of the distance modulus, before 
folding in an explicit prior on VAi, but after accounting for the priors on M and [Fe/H] 
(Eq. 3). In this section we present some example =Sf (PA^), but foremost show what the 
systematic impact on is of neglecting the M and [Fe/H] priors. 

When estimating the distance to a given star, without the benefit of external prior 
information, one would evaluate Eq. 3 presuming that Pprior(M) and Pprior([Fe/H]) are 
constant. 

To test the impact of Ppi.ior(M), we estimate the distances for two cases. No priors 
applied (pprior(M) = 1 and Pprior([Fe/H]) = 1); and only the prior on the luminosity function 
applied (p(L) oc L~^-^ and Ppnor([Fe/H]) = 1). The distance modulus estimate that neglects 
the explicit priors is denoted as DMq, while the distance modulus with only Pprior(M) applied 
is marked as DMl. The top panel of Figure M illustrates the importance of including the 
'luminosity prior', by showing the systematic difference in VAi that results from neglecting 
it. For stars near the tip of the giant branch the bias is very small, but for stars near the 
bottom of the giant branch, the mean systematic bias of neglecting the prior information is 
0.08 mag, with systematic bias of ~ 0.2 mag in some cases. 
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To test the impact of Pprior([Fe/H]), we estimate the distances where only the metaUicity 
prior was apphed, and mark the relevant VAi as DM[fc/h]- Compared with the distance 
modulus with no priors applied, DMq, we find Ppi.ior([Fe/H]) can correct a mean overestimate 
of 0.03 mag on the VA4 for the metal-poor stars and a mean underestimate of 0.02 mag on 
the T>Ai for the metal-rich ones, but the neglect of the metaUicity prior causes a smaller 
bias in VAi than neglecting the luminosity prior (0.05 vs. 0.2 mag at maximum), as shown 
in Figure [9] (middle panel). The distance modulus bias caused by neglecting both priors is 
presented in the bottom panel of Figure [9l Neglecting both priors causes a mean bias of 0.08 
mag and a maximum bias of 0.2 mag in distance modulus. 

4.3. Distance Precision Test using Clusters 

We use known members from three clusters (M92, M13 and NGC6791) to test the 
estimates of the VAi, because they have more spectroscopic member giants than M71. 
Figure [TD] shows all the member stars in each cluster and the member giants used for the 
test are highlighted. The Bayesian method was applied to estimate VAi for each member 
giant, adopting the same luminosity prior, but a Gaussian distribution of [Fe/H] centered 
at the literature cluster metaUicity, and with the same [Fe/H] dispersion as known cluster 
members, as the metaUicity prior. Figure [TT] shows our estimates of VAi are consistent with 
the literature VAicc (shown in Table 1), derived by main-sequence fitting, to within 1 a. 

5. Summary and Conclusions 

We have implemented a probabilistic approach to estimate the distances for SEGUE K 
giants in the Galactic halo. This approach folds all available observational information into 
the calculation, and incorporates external information through priors, resulting in a VAi 
likelihood for each star that provides both a distance estimate and its uncertainty. 

The priors adopted in this work are the giant-branch luminosity function derived from 
globular clusters, and the ensemble metaUicity distributions for different SEGUE K-giant 
target categories. We show that these priors are needed to prevent systematic overestimate 
of the distance moduli by up to 0.2 mag. The role of the priors are important, and make 
the results more reasonable. 

We use empirical color-magnitude fiducials from old stellar clusters to obtain the pre- 
dicted colors c(M, [Fe/H]), which are needed to calculate ^(VAi). Ultimately, the best 
estimates of the distance moduli and their errors can be estimated using the peak and central 
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68% interval oi^{VM). 

With this approach we obtain Bayesian estimates of distance moduh, and thus, absolute 
magnitudes and distances, for 4781 K giants that have a mean distance precision of 12%, 
or ±0.25 mag in VAi and M^. The sample contains 270 stars beyond rcc = 50 kpc, which 
makes it by far the largest sample of distant stellar halo stars with measured radial velocities 
and distances to date. 

We present an online catalog containing the distance moduli, observed information 
and SSPP atmospheric parameters for the 4781 SEGUE K giants. For each object in the 
catalog we also list some of the basic observables such as (RA,Dec), extinction-corrected 
apparent magnitudes and colors, as well as the information obtained from the spectra- 
heliocentric radial velocities plus SSPP atmospheric parameters. In addition, we provide 
the Bayesian estimates of the distance moduli, distances to the Sun, Galactocentric dis- 
tances, the absolute magnitudes and their uncertainties, along with the distance moduli at 
(5%, 16%, 50%, 84%, 95%) confidence of ^{VM) . 

We tested the accuracy of our distance estimates using giants in the old clusters M13, 
M92 and NGC6791, which have SEGUE observations. We found that the distance estimates 
for the individual cluster member stars derived with our approach are consistent with the 
literature distance moduli of the clusters. 

We caution the reader that the n(d, [Fe/H]) cannot be used to obtain the halo profile and 
the metallicity distribution directly, because the complex SEGUE selection function needs 
to be taken into account. 
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Fig. 1. — : Color-color diagram for confirmed K giants in our sample, with DR9 SPP estimates 
for [Fe/H] color coded as shown in the vertical bar. 
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Fig. 2. — : Distribution of our K-giant sample in the color- magnitude and color- metallicity 
plane. It can be seen that the most common stars are the intrinsically fainter blue giants, 
as we would expect from the giant-branch luminosity function. The possible HB/RC stars 
are over-plotted as black dots, and the relation between [Fe/H] and {g — r)^^ is marked as 
the black solid line. 
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M 

Fig. 3. — : Luminosity functions for the giant branch of two globular clusters in differ- 
ent bands, compared with two theoretical giant-branch luminosity functions. We scaled 
logio(dN) to make the luminosity functions separate from each other, because the slope 
is the important parameter to test whether the luminosity functions are consistent. All 
the luminosity functions follow power laws with nearly the same slope of 0.32, so that the 
theoretical and observational luminosity functions are consistent, and both are insensitive 
to changes in metallicity and passbands. Therefore, the prior probability adopted for the 
absolute magnitude in the analysis is p(M) = 10'^'^^^/17.788, which is based on a vari- 
ety of theoretical and empirical giant-branch luminosity functions, and whose integral over 
[-3. 5, +3. 5] has been normalized to unity, 'k' stands for the slope. 
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Fig. 4. — : (Upper panel) Variation of the metallicity distribution with apparent tq magnitude 
for four sub-categories. (Lower panel) The four [Fe/H] priors adopted in this work. The 
integral of p([Fe/H]) has been normalized to unity. 
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Fig. 5. — : Interpolation of the four giant-branch fiducials (thick hnes), to obtain the relation, 
r)o = f{Mr, [Fe/H]), for any set of (M^, [Fe/H]). The thin lines show a set of interpolated 
fiducials. No values outside the extreme fiducials are used. 
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Fig. 6. — : Examples of =Sf(PA^) for three stars, with and without accounting for the 
metaUicity and luminosity- function priors (see §4.2). The black line indicates the most likely 
VAi under the three assumptions, showing the systematic overestimate of the V/A if priors 
are neglected. 
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|\/|^ logio(N) 

Fig. 7. — : Results of the distance estimates for 4781 K giants. (Upper panel) The distribution 
of the relative error in distance vs. distance. (Middle panel) The distribution of the error 
in absolute magnitude vs. absolute magnitude. (Lower panel) The distribution of the error 
in distance moduli vs. distance moduli. Note that the fractional distance estimates are less 
precise for nearby stars, because the lower part of the giant branch (less luminous, therefore 
more nearby) is steep in the color-magnitude diagram, particularly at low metallicities (see 
Figure El) . 
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Fig. 8. — : (Upper panel) The distribution of K giants on the CMD plot. (Lower panel) 
The distribution of the mean error in the absolute magnitude, M^, as a function of and 
((? — r)o. The upper panel shows that the sample contains a large fraction of relatively nearby 
giants of moderate luminosity (M^ ~ 0). (Lower panel) The luminosity estimates for stars 
in lower part of the CMD are less precise, because the isochrones are much steeper in this 
part, especially for low metallicities. 
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Fig. 9. — : The distance modulus bias caused by neglecting the priors on the luminosity 
function and metallicity distribution of the K giants. The luminosity function prior can help 
correct a mean overestimate of 0.1 mag in the distance moduli, and a maximum overestimate 
of ~ 0.2 mag in some cases. While the [Fe/H] prior's impact is much smaller, it can help 
correct a mean of 0.02 mag overestimate or underestimate on the distances in the metal-poor 
or metal-rich tails. The bottom panel shows the total impacts of the luminosity prior and 
the metallicity prior. 
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Fig. 10. — : The color-magnitude diagrams for the three clusters used for distance precision 
test. The black dots are cluster member stars that we used to derive the fiducials, and the 
red dots are the K giants we used to examine the Bayesian approach. 
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Fig. 11. — : The difference between our estimate of the Clusters' VAi and the hterature 
VAi vs. ((yf — r)o. Black triangles indicate the case using the cluster metallicity priors, which 
indicates that the Bayesian estimates of VAi are consistent with the literature VAicc, to 
within 1 a. 
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Table 1:: Parameters of the Fiducial Clusters 



NGC Messier E(B-V) (m - M)o [Fe/H] 



6341 


M92 


0.02^ 


14.64^ 


-2.38^ 


6205 


M13 


0.02^ 


14.38'^ 


-1.60^ 


6838 


M71 


0.28= 


12.86'= 


-0.81^ 


6791 




0.16<i 


13.01<i 


+0.39'^ 



" IKraft fc Ivana (j2003l ): their globular cluster metallicity scale is based on the Fell lines from high-resolution 
spectra of giants. 

' ICarretta et"all |2000l ): (m — M)q derived from the Hipparcos sub-dwarf fitting. 

iGrundahl et al.l (|2002l ): (m — A/)o derived from the Hipparcos (jPerrvman et al.lll997l ) sub-dwarf fitting. 
iBrogaard et al.l (j201l[ ): (m — M)o is based on (m — M)„ assuming = 3.1E{B — V). 



simple average of the [Fe /H] + 0.29, +0A7, +0A an d +0.39 bv lBrogaard et al.l (|201ll) . iGratton et al.l (|2006f ). 
Peterson fc GreenI |l998h . and lCarraro et all |2006l) respectively. 



Table 2:: Metallicity and Color of the Red Hori z ontal- 



Branch Onset for the Eight Clusters in I An et al.l (120081 ) 



Name of Clusters [Fe/H] {g - r)^^ 

NGC6791 +0.39 1.13 

M5 -1.26 0.61 

M3 -1.50 0.59 

M13 -1.60 0.58 

M53 -1.99 0.54 

M92 -2.38 0.53 

M15 -2.42 0.53 

The first column are the names of the clusters and 
the next two columns contains the [Fe/H] of the clus- 
ters and the extinction corrected color [g — r)o. 
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Table 3:: Interpolated Fiducial at 



[Fe/H]=-1.18 



(g - r)o 




[Fe/H] 


0.47 


3.00 


-1.18 


0.52 


2.01 


-1.18 


0.57 


1.31 


-1.18 


0.61 


0.68 


-1.18 


0.66 


0.14 


-1.18 


0.70 


-0.30 


-1.18 


0.75 


-0.65 


-1.18 


0.79 


-0.94 


-1.18 


0.84 


-1.18 


-1.18 


0.88 


-1.39 


-1.18 


0.93 


-1.58 


-1.18 


0.98 


-1.75 


-1.18 


1.02 


-1.90 


-1.18 


1.07 


-2.04 


-1.18 


1.11 


-2.19 


-1.18 


1.16 


2.33 


-1.18 


1.20 


-2.45 


-1.18 


1.25 


-2.55 


-1.18 


1.29 


-2.65 


-1.18 


1.34 


-2.74 


-1.18 



This is an example of one interpo- 
lated fiducial with [Fe/H] =-1.18. 
The entire catalog of 20 interpo- 
lated fiducials with [Fe/H] rang- 
ing from -2.38 to +0.39 is pro- 
vided in the electronic edition of 
the Journal. 



Table 4:: List of 4781 K Giants Selected from SDSS DR9 



RA(J2000) Dec(J2000) rg Aro (9 - r)o A(g - r)o RV ARV Tett [Pc/H] log g BA^poak 

(dcg) (dcg) (mag) (mag) (mag) (mag) (kma^-^) (kms^'-) (K) (mag) (mag) 



359.5853 36.4018 15.139 0.010 0.713 0.022 -71.9 2.1 4979 -2.17 1.80 16.59 16.15 

358.5352 36.2383 16.762 0.010 0.623 0.014 -198.4 3.5 5140 -1.32 2.40 16.42 16.08 

358.2261 36.1001 16.002 0.020 0.923 0.028 -60.4 4.5 4496 -2.78 1.10 18.59 18.39 

357.3861 35.7609 16.239 0.010 0.647 0.022 -172.6 3.1 5062 -1.48 2.00 16.33 15.84 

357.7179 37.0025 15.519 0.010 0.685 0.014 -63.4 2.3 4936 -2.10 1.50 16.71 16.33 



RA(J2000) 


Doc(J2000) 










AT>M 


Mr 


AJ\/r 


d 


Ad 


I'GC 


Aroc 


(deg) 


(dcg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


359.5853 


36.4018 


16.32 


16.54 


16.73 


16.84 


0.21 


-1.45 


0.21 


20.79 


1.94 


24.59 


1.84 


358.5352 


36.2383 


16.21 


16.41 


16.61 


16.73 


0.20 


0.34 


0.20 


19.22 


1.75 


23.01 


1.65 


358.2261 


36.1001 


18.47 


18.58 


18.68 


18.73 


0.10 


-2.59 


0.10 


52.35 


2.41 


55.34 


2.39 


357.3861 


35.7609 


16.03 


16.31 


16.57 


16.73 


0.27 


-0.09 


0.27 


18.48 


2.25 


22.20 


2.11 


357.7179 


37.0025 


16.48 


16.68 


16.84 


16.94 


0.18 


-1.19 


0.18 


22.01 


1.80 


25.62 


1.72 


The first 


two columns contains the positi 


ion (RA, Dec) for each object. 


The mag 


nitude, color, and 


their errors are 


in the next 


four columns: 


corrected for 


extinction. The 


I>JM at the 



peak and (5%, 16%, 50%, 84%, 95%) confidence of ^ (X^A^) are listed in the next six columns. T^-Mpeak ^^e best estimate of the distance modulus for the K giant. The ^q^j^ is the 
uncertainty of the distance modulus, which is calculated from {T^-M.q4% — '^-'^ 16% )/-^- The last six columns are absolute magnitude and distances calculated from 'DA4pgg|_k^ assuming 

■DA4 + 5 

Rq — 8.0 kpc {i.e. M,- — Tfl — A4pcak 1 d — 10 5 ). The complete version of this table is provided in the electronic edition of the Journal. The printed edition contains only a 
sample. 



